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Abstract 

The charm fragmentation function has been measured in D* photoproduction 
with the ZEUS detector at HERA using an integrated luminosity of 120 pb^^. 
The fragmentation function is measured versus z = {E + p\\)^'/20'=\ where 
E is the energy of the D* meson and p\\ is the longitudinal momentum of the 
D* meson relative to the axis of the associated jet of energy i?-""*. Jets were 
reconstructed using the clustering algorithm and required to have transverse 
energy larger than 9 GeV. The D* meson associated with the jet was required 
to have a transverse momentum larger than 2 GeV. The measured function is 
compared to different fragmentation models incorporated in leading-logarithm 
Monte Carlo simulations and in a next-to-leading-order QCD calculation. The 
free parameters in each fragmentation model are fitted to the data. The extracted 
parameters and the function itself are compared to measurements from e"'"e~ 
experiments. 
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1 Introduction 



The production of a charm hadron is described as the convolution of the perturbative 
production of a charm quark and the non-perturbative transition of a charm quark to 
a hadron. The non-perturbative component is assumed to be universal, i.e. indepen- 
dent of the initial conditions. It is described by so-called fragmentation functions which 
parametrise the transfer of the quark's energy to a given hadron. The free parameters 
are determined from fits to data. The transition of a charm quark to a D* meson is the 
subject of this paper. 

The parameters of the various fragmentation function ansatze were so far derived from 
data obtained at e^e~ colliders. The e^e~ data span a wide range of centre-of-mass 
energies and the fragmentation of a charm quark to a D* meson has been measured 
many times [1], most recently by the CLEO [2] and Belle [3] collaborations at a centre- 
of-mass energy of ~10.5GeV and the ALEPH [4] collaboration at 91.2 GeV. Due to 
scaling violations in QCD, the dependence of the fragmentation function on production 
energy [1,5] is expected to follow the DGLAP equations [6]. 

The fragmentation function has recently been measured by the HI Collaboration for the 
production of D* mesons in deep inelastic scattering (DIS) [7]. A measurement of the 
fragmentation function at HERA and its comparison with that deduced from experi- 
ments at e^e~ colliders provides a measure of the universality of charm fragmentation 
and further constrains its form. The analysis presented here has been performed in the 
photoproduction regime in which a quasi- real photon of low virtuality, Q^, is emitted from 
the incoming electron or positron and collides with a parton in the proton. 



2 Experimental conditions 

The analysis was performed using data collected with the ZEUS detector at HERA 
during 1996-2000. In this period, HERA collided electrons or positrons with energy 
= 27.5 GeV and protons with energy Ep = 820 GeV (1996-1997) or Ep = 920 GeV 
(1998-2000) corresponding to integrated luminosities of 38.6 ± 0.6 and 81.9 ± 1.8 pb~^ 
and to centre-of-mass energies y/s = 300 GeV and y/s = 318 GeV, respectively. 

A detailed description of the ZEUS detector can be found elsewhere [8]. A brief outline 
of the components that are most relevant for this analysis is given below. 

Charged particles were tracked in the central tracking detector (CTD) [9], which operated 
in a magnetic field of 1.43 T provided by a thin superconducting coil. The CTD consisted 
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of 72 cylindrical drift chamber layers, organised in 9 superlayers covering the polar-angle^ 
region 15° < 9 < 164°. The transverse-momentum resolution for full-length tracks was 
(^{Pt)/pt = 0.0058pT © 0.0065 © 0.0014/pt, with pr in GeV. 

The high-resolution uranium-scintillator calorimeter (CAL) [10] consisted of three parts: 
the forward (FCAL), the barrel (BCAL) and the rear (RCAL) calorimeters. Each part was 
subdivided transversely into towers and longitudinally into one electromagnetic section 
(EMC) and either one (in RCAL) or two (in BCAL and FCAL) hadronic sections (HAC). 
The smallest subdivision of the calorimeter was called a cell. The CAL energy resolutions, 
as measured under test-beam conditions, were a{E)/E = 0.18/V^ for electrons and 
a{E)/E = 0.35/ VE for hadrons, with E in GeV. 

The luminosity was measured from the rate of the bremsstrahlung process ep —>■ e'-fp, 
where the photon was measured in a lead-scintillator calorimeter [11] placed in the HERA 
tunnel at Z = — 107 m. 



3 Event selection and reconstruction 

A three-level trigger system was used to select events online [8,12,13]. At the first- and 
second-level triggers, general characteristics of photoproduction events were required and 
background due to beam-gas interactions rejected. At the third level, a version of the 
tracking information close to the offline version was used to select D* candidates. 

Kinematic variables and jets were reconstructed offline using a combination of track and 
calorimeter information that optimises the resolution of reconstructed kinematic vari- 
ables [14]. A selected track or calorimeter cluster is referred to as an Energy Flow Object 
(EFO). The jets were reconstructed with the /c^ cluster algorithm [15] in its longitudi- 
nally invariant inclusive mode [16], where the parameter R is chosen equal to 1. Jets 
were formed from the EFOs with at least one jet required to have transverse energy, 
E^r^^ > 9 GeV and pseudorapidity, {rj^'^^l < 2.4. The photon-proton centre-of-mass energy, 
Wyp, was calculated using the formula Wyp = ^J'^Epi^i ~ Pz,i), where the sum runs 
over the energy and longitudinal momentum component of all EFOs. Due to trigger re- 
quirements and beam-gas background at low W^p and background from DIS events at 
high W-yp, the requirement 130 < W^p < 280 GeV was made. Neutral current DIS events 



The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction" , and the X axis pointing left towards 
the centre of HERA. The coordinate origin is at the nominal interaction point. The pseudorapidity 
is defined as i] = — In (tan |), where the polar angle, 6, is measured with respect to the proton beam 
direction. 
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with a scattered electron or positron candidate in the CAL were also removed by cut- 
ting [17] on the inelasticity, which is estimated from the energy, E'^, and polar angle, 
6'g, of the scattered electron or positron candidate using 1/^ = 1 — — cos6'g). Events 
were rejected if ?/e < 0.7. 

The D* mesons were identified using the decay channel D*^ D^T^t with the subsequent 
decay K^n^ and the corresponding anti-particle decay. They were reconstructed 

from charged tracks in the CTD using the mass-difference technique [18]. Tracks with 
opposite charges and transverse momenta greater than 0.5 GeV were combined into pairs 
to form D° candidates. No particle identification was used, so kaon and pion masses were 
assumed in turn for each track to calculate the invariant mass M^K-k). A third track, 
assumed to be the soft pion, 7r+, with transverse momentum greater than 0.12 GeV and 
of opposite charge to the kaon, was combined to form a D* candidate with invariant mass 
M^KiTTis). The D* candidates were then required to have > 2 GeV and \ri^ \ < 1.5. 

To minimise background, narrow windows were selected for the mass difference, AM = 
M{K'K'Ks) - M{KTr), and the mass of the D° meson: 0.1435 < AM < 0.1475 GeV and 
1.83 < M{D^) < 1.90 GeV. For background determination, candidates with wrong- 
charge combinations, in which both tracks forming the candidates have the same 
charge and the third track has the opposite charge, were also retained. The same kinematic 
restrictions were applied as for those candidates with correct-charge combinations. The 
normalisation factor of the wrong-charge sample (a value of 1.02 for the distribution after 
all requirements shown in Fig. 1) was determined as the ratio of events with correct-charge 
combinations to wrong-charge combinations in the region 0.150 < AAf < 0.165 GeV. A 
cut of pJ^* I E^^^^° > 0.1 was imposed to further reduce combinatorial background, where 
^e^>io° ^Yi^ transverse energy measured using all EFOs outside a cone of 10° in the 
forward direction. The forward region was excluded because of the strong influence of the 
proton remnant [19]. 

Finally, the D* meson was associated with the closest jet (with E^^^ > 9 GeV and 
jiyjetj ^ 2.4) in ?7 — space and requiring R ^= a/ {rj^'^^ — rj^'Y + {(j)^'^^ — (f)^*Y^ < 0.6. 

The combined efficiency for all the above requirements was about 35%. A clear D* mass 
peak above a relatively small background is shown in Fig. 1. Subtraction of the back- 
ground of 634 ±30 candidates, estimated from the wrong-charge sample, gave 1307 ±53 
D* mesons. The background was subtracted bin-by-bin as a function of the measured 
fragmentation variable and all other subsequent distributions. 
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4 Fragmentation variables and kinematic region 



In e~^e~ collisions, at leading order (LO), the two produced charm quarks each carry half 
of the available centre-of-mass energy, ^/s. The fragmentation variable of a D* meson can 
therefore be simply related to one of the two produced jets. In ep collisions, the definition 
of the fragmentation variable is not so simple as only a fraction of the available centre- 
of-mass energy contributes to the production of charm quarks in the hard scattering 
process. However, charm quarks produced in the hard scatter form final-state jets of 
which the meson is a constituent. Therefore, the fragmentation variable, z, is calculated 
as z = {E + p\\)^ /{E + P\\)^'^^i where p\\ is the longitudinal momentum of the D* meson 
or of the jet relative to the axis of the associated jet of energy, E^^^, where all quantities 
are given in the laboratory frame. As the jets are reconstructed as massless objects, z 
simplifies to: 



The analysis was performed in the photoproduction regime with 130 < W^p < 280 GeV 
and < 1 GeV^. The D* meson was required to be in the region |?7^*| < 1.5 and 
pJ^* > 2 GeV. The D* meson was included in the jet-finding procedure and was thereby 
uniquely associated with one jet only. Each jet associated with a D* was required to 
satisfy < 2.4 and E^^' > 9 GeV. 

Cuts on the minimum jet transverse energy and minimum D* transverse momentum 
will lead to a bias in the z distribution as 2; ~ /E^°^. Therefore the minimum jet 
transverse energy was chosen to be as high as possible and the minimum D* transverse 
momentum to be as low as possible whilst maintaining statistical precision. With the 
above requirements, the z distribution is unbiased above 0.22. 

5 Fragmentation models 

Various parametrisations of fragmentation functions have been proposed. Those consid- 
ered in this paper are detailed below. 

A parametrisation often used to describe the fragmentation of heavy quarks is the function 
from Peterson et al. [20] which has the form 




D* 
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(2) 




where e is a free parameter. 
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The function from Kartvelishvili et al. [21] has the form 



(3) 



where a is a free parameter. 

Within the framework of the Artru-Mennessier model [22], Bowler [23] developed a frag- 
mentation function for heavy quarks of mass, tuq, which has the form 



where a and b are free parameters and tq is predicted to be unity. The quantity m± is the 
transverse mass of the hadron, = im? + (p^')^, where m is the hadron's mass and p^^^ 
the transverse momentum relative to the direction of the quark. The additional freedom 
given by rg allows a smooth transition to the symmetric Lund form [24] [rq = 0) used to 
describe light-quark fragmentation. 

6 Monte Carlo models 

Monte Carlo (MC) models were used both to calculate the acceptance and effects of the 
detector response and to extract fragmentation parameters. The programmes Herwig 
6.1 [25] and Pythia 6.1 [26] which implement LO matrix elements followed by parton 
showers and hadronisation were used to model the final state. Different parameter settings 
were used in the MC models when correcting the data or when extracting fragmentation 
parameters; the settings used when extracting fragmentation parameters are given in 
Section 8.1. The MC used to correct the data had the default settings, apart from the 
following changes: the fraction of charged D mesons produced in a vector state was set 
to 0.6 [27]; and the excited D-meson production rates were set to non-zero values^ [28]. 

The ZEUS detector response was simulated in detail using a programme based on Geant 
3.13 [29]. The Pythia 6.1 MC programme was used with two different fragmentation 
schemes: the default which is the Lund string model [30] modified according to Bowler 
for heavy quarks; and the Peterson fragmentation function with e = 0.06 (see Section 5). 
The Herwig 6.1 MC programme uses a cluster model [31] for its fragmentation. 



^ These changes correspond to the Pythia parameters: PARJ(13) = 0.6, PARJ(14) = . 13, PARJ (15) 
= 0.01, PARJ(16) = 0.03 and PARJ (17) = 0.13. 
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7 Data correction and systematic uncertainties 



The data were corrected for acceptance and effects of detector response using a bin-by- 
bin metliod with the Pythia simulation used as the central MC. The distribution of the 
difference in z between hadron and detector levels is symmetric, has a mean of zero and 
a width of 0.06 units. To maintain both high purity and the validity of the bin-by-bin 
method, a bin width of at least twice this value (0.14 units) was chosen. The rate due 
to h quarks was subtracted using the Pythia MC prediction normalised to a previous 
measurement of jet photoproduction [32]. Therefore the cross section as a function of z is 
for processes in which an initial-state charm quark hadronises to a D* meson. A detailed 
analysis [33] of the possible sources of systematic uncertainty was performed. The sources 
are: 

b\ the use of an alternative fragmentation model in the Pythia MC simulation (see 
Section 6). As the Herwig MC simulation gave a poor description of the data, it was 
not used to correct the data; 

^2 the h fraction subtracted was changed by (a) -|-30% and (b) —30% in accordance 
with the level of agreement between data and Pythia MC predictions [32] for jet 
photoproduction; 

53 the total energy in the jet reconstructed from the CAL EFOs was varied by (a) +3% 
and (b) —3% in the simulation, in accordance with the uncertainty in the jet energy 
scale; 

^4 the range of W^^ was changed to (a) 124 < W^^ < 267 GeV and (b) 136 < W^p < 
293 GeV, in accordance with the resolution; 

^5 the cut on i?;^* was changed to (a) 10 GeV and (b) 8 GeV, in accordance with the 
resolution; 

^6 the value of the cut on p!j^* / E^^^° was varied to (a) 0.08 and (b) 0.12; 

^7 the lower (upper) bound on the normalisation region for the wrong-charge candidates 
was changed to (a) 0.152 ((b) 0.163) GeV. 

The cuts on 77-^°*, 77^ and pf^ were also varied in accordance with their resolution and 
produced negligible effects. The values of the above uncertainties for each bin in the 
normalised cross section, {l/a)da/dz, are given in Table 1. The individual systematic 
uncertainties were added in quadrature separately for the positive and negative deviations 
from the nominal cross-section values to obtain the total systematic uncertainties. The 
systematic uncertainties on the fits of the various fragmentation parametrisations to the 
data described in Section 8 were obtained from fits to the cross section for each systematic 
variation. The resulting variations in a given fragmentation parameter were added in 
quadrature to yield the systematic uncertainty on that parameter. 
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8 Results 



The distributions of the variables 2, , where is the transverse momentum of the 
D* meson relative to the jet, p^*, r]^* , E}^^ and "rf^^ are shown in Fig. 2 and compared to 
the distributions from the MC programmes, normalised to the data. Also shown is the 
prediction of the Pythia simulation for the production of beauty quarks subsequently 
producing a D* meson; this amounts to about 6%. The z distribution is reasonably well 
described by the Pythia MC predictions, whereas the Herwig prediction does not de- 
scribe the data. This can be seen in the differences between the measured p^* distribution 
and that predicted by Herwig. The MC predictions for the E^^^ distribution are, how- 
ever, similar and agree reasonably well with the measurement. For the p^^' distribution, 
the Pythia simulations give a good description of the data and are again better than 
that from Herwig. This shows that the Pythia MC model using both the Bowler and 
Peterson fragmentation for charm quarks gives a good description of the transverse as well 
as the longitudinal component of the D* fragmentation process. The distribution of the 
pseudorapidities of both jet and D* are similarly well described by both MC programmes. 
As the Herwig MC model is known to give a better description than Pythia of data [34] 
sensitive to the parton-shower model, the differences shown here suggest that the cluster 
model does not describe the hadronisation process of charm quarks to D* mesons. 

The normalised differential cross section, \j aida j dz\ is presented in the kinematic region 
< 1 GeV^ and 130 < W^^ < 280 GeV, requiring at least one jet with E^r^' > 9 GeV and 
|?7j'^*| < 2.4. A D* meson with p^* > 2 GeV and \ri^* \ < 1.5 was required to be associated 
with any jet that satisfied the above jet requirements on E^r^^ and rj^^^. The D* meson 
was included in the jet-finding procedure and was thereby uniquely associated with one 
jet only. The values of the cross section are given in Table 2 and shown in Figs. 3 and 4 
compared to various expectations. In Fig. 5, the same data are shown compared with 
results from e~^e~ experiments. 

8.1 Comparison with fragmentation models in Pythia 

The normalised cross section is shown in Fig. 3 compared to the Pythia MC simulation 
using different fragmentation models. The original default settings for Pythia 6.1 were 
used with the proton and photon parton density functions set to GRV94 LO [35] and GRV- 
LO [36], respectively and a different value for the maximum parton virtuality allowed in 
space-like showers (PARP(67) in Pythia changed from 1.0 to 4.0 [37]). Otherwise, only 
the fragmentation parameters considered (see Section 5) were varied. 

The default fragmentation setting in the simulation is the symmetric Lund string frag- 
mentation modified for heavy quarks according to Bowler (see Eq. 4). Three predictions 
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for different values of rq are shown compared to the data in Fig. 3(a). The default predic- 
tion with tq = 1 gives a reasonable description of the data; as tq decreases, the prediction 
deviates more and more from the data. 

The Peterson function (see Eq. 2) and the option to vary e is available within the Pythia 
simulation. The value of e was varied in the range 0.01 to 0.1, with the Lund string 
fragmentation model used for lighter flavours. For each value in the MC simulation, the 
full event record was generated and the kinematic requirements applied, allowing a direct 
comparison to the data. The result of varying e is shown in Fig. 3(b). Here it can be seen 
that values as low as e = 0.01 are disfavoured, producing a much harder spectrum than 
the data, while values as high as e = 0.1 result in too soft a spectrum and are therefore 
also disfavoured. The result of fitting the MC to the data was e = 0.062 ± 0.007l[]:gg^ 
where the first uncertainty is statistical and the second systematic. The value is consistent 
with the default value in the MC of e = 0.05 which was obtained from comparisons [26] 
with LEP and SLD data at the mass. The fitted value was then used in the MC and 
the result compared in Fig. 3(b); the data are well described. 



8.2 Comparison with next-to- leading-order QCD calculations 

The data were compared with a next-to-leading-order (NLO) QCD prediction [38] which 
is a fixed-order calculation from Frixione et al. (FMNR). As default, the programme is 
interfaced to the Peterson fragmentation function; the function from Kartvelishvili et 
al. (see Eq. 3) was also implemented. The other parameters used in the NLO QCD 
calculation were as follows: the renormalisation and factorisation scales were set to /i = 
■\/ {{PtY) + ml, where {{PtY) is the average squared transverse momentum of the two 
charm quarks and rric = 1.5 GeV; the proton parton density function was CTEQ5M1 [39]; 
and the photon parton density function was AFG-HO [40]. 

As the final state particles in the NLO QCD calculation are partons, to enable a fair 
comparison with the data, the predictions were corrected for effects of hadronisation 
using a bin- by-bin procedure according to Act = Acx'^^^ ■ Chad? where Aa^^'~' is the 
cross section for partons in the final state of the NLO calculation. The hadronisation 
correction factor, Chad? was defined as the ratio of the cross sections after and before the 
hadronisation process, Chad = Af^uc^"^^^ / '^^mc°^^ ^ where the partons used are those after 
parton showering. The values of Chad from Pythia were used for the central results. As 
the results of Herwig do not describe the data (see Section 6), they are used only as a 
systematic check. The prediction from this combination of NLO QCD and hadronisation 
correction is termed "FMNRxChad"""" • The values of Chad are given for Pythia and 
Herwig in Table 2. 
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The result of varying e in the Peterson function and a in the Kartvehshvih function for the 
predictions of FMNRxC^g^""" ^^e shown in Figs. 4(a) and (b), respectively. The data again 
show sensitivity to these fragmentation functions and can constrain their free parameters. 
The results of fits to the data are e = 0.079 ± 0.008l|J:[JJ^ and a = 2.67 ± O.lStoH for the 
Peterson and Kartvelishvili functions, respectively, where the first uncertainty is statistical 
and the second systematic. 

A number of parameter settings which are commonly used in comparison with data [34] 
were considered. Using Chad from Herwig gave e = 0.094 ± 0.008 and a = 2.46 ± 0.17, 
where the uncertainty is statistical only. The effect of the input parameters in the NLO 
QCD programme was checked by changing the renormalisation scale and charm mass 
simultaneously to 2/i and 1.7 GeV and 0.5/i and 1.3 GeV. The different settings gave 
values of e [a) of 0.082 (2.55) and 0.077 (2.80), respectively; the uncertainty from the 
NLO QCD input parameters is significantly smaller than the experimental uncertainties. 

The default e value used so far in NLO QCD calculations, extracted from a fit [41] to 
ARGUS [42] data, was 0.035. As the perturbative part of the production in calculations 
of e~^e~ and ep cross sections depends on the scale of the process and colour connections 
between the outgoing quarks and the proton remnant can have an effect, the values of e 
extracted with NLO QCD from e~^e~ and ep data may not necessarily be the same. This 
illustrates that care is needed in choosing the appropriate fragmentation parameter. 

8.3 Measurement of (z) and comparisons with e+e~ data 

In Fig. 5, the ZEUS data are shown compared with measurements from the Belle [3], 
CLEO [2] and ALEPH [4] collaborations in e^e~ interactions. The Belle and CLEO 
data are measured at a similar centre-of-mass energy of about 10.5 GeV, whereas the 
ALEPH data was taken at 91.2 GeV. The corresponding scale of the ZEUS data is given 
by twice the average transverse energy of the jet, 23.6 GeV, and is between the two e'^e~ 
centre-of-mass energies. 

Although using a different definition for z, the general features of the data presented here 
are similar to those at e^e~ experiments. However the ZEUS data are shifted somewhat to 
lower values of z compared to the CLEO and Belle data with the ALEPH data even lower. 
This can be seen more quantitatively by extracting the mean value of the distribution, 
(z) = 0.588 ± 0.025 (stat.) ± 0.029 (syst.). The Pythia MC programme was used to 
extrapolate the phase space to p^* = and to correct for the subsequent exclusion of 
the region 0<2;<0.16. It was also used to correct for the finite bin size. The resulting 
factor was 0.961. The corrected value. 



9 



{z) = 0.565 ± 0.024 (stat.) ± 0.028 (syst.) 



(5) 



and those from ALEPH, Belle and CLEO are shown in Table 3. It should be noted that 
the ALEPH data uses the beam energy as the scale rather than the jet energy which, 
due to hard gluon emission, would be a better equivalent to the jet energy used in this 
analysis. The usage of jet energy for ALEPH data would lead to an increase in (z). 
Although the uncertainties on the current measurement are larger than those from the 
e"'"e~ experiments, the value is qualitatively consistent with expectations from scaling 
violations in QCD in which {z) decreases with increasing energy [43]. 

9 Summary 

The fragmentation function for D* mesons has been measured in photoproduction at 
HERA using the variable z = {E + p\\)^* /2E^^^ and requiring a jet with E^^* > 9 GeV and 
I'tf'^^l < 2.4 to be associated with a D* meson in the range p^* > 2 GeV and |77^*| < 1.5. 

The data are compared to different fragmentation models in MC simulations and a NLO 
QCD calculation. The cluster model used in the Herwig programme does not describe 
the data. Within the framework of NLO QCD and the Pythia simulation, the free 
parameters of the Peterson fragmentation function and, for NLO QCD, the Kartvelishvili 
function have been fitted. 

The value of e in the Peterson function, extracted within the framework of NLO QCD, is 
different to that extracted using data from e^e~ collisions. As the perturbative aspects 
of the corresponding calculations and the energy scales are different, the results are not 
expected to be the same. Future calculations of charm hadron cross sections at NLO QCD 
at HERA should always use the appropriate values. Within the consistent framework 
given by the Pythia model, the extracted fragmentation parameters agree with those 
determined in e~^e~ data. 

The fragmentation function and the (z) are different to those measured at different centre- 
of-mass energies in e+e~ collisions; the measured {z) is higher than the ALEPH data and 
lower than the CLEO and Belle data, qualitatively consistent with the scaling of this 
variable as predicted by QCD. 
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Table 1: Individual sources of systematic uncertainty (in %) per bin of the nor- 
malised cross-section {\/(j)da/dz. The description of each variation is given in 
Section 7. 
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Table 2: Measured normalised cross-section {1/ a)da / dz. The statistical ((^statj CLnd 
systematic (Ssyst) uncertainties are shown separately. The bin-by-bin corrections for 
hadronisation (see Section 8.2) are shown for Pytri A, C^g^^ , and Herwig, C™^. 
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Collaboration 


Scale (GeV) 


Measured variable 


{z) ± stat. ± syst. 


ALEPH 


91.2 




0.4878 ±0.0046 ±0.0061 


Belle 


10.6 




0.61217 ±0.00036 ±0.00143 


CLEO 


10.5 




0.611 ±0.007 ±0.004 


ZEUS 


23.6 


((E+P||)^72Ej^*) 


0.565 ±0.024 ±0.028 



Table 3: Mean value, (z) , of the fragmentation function in e~^e~ collisions, 
ALEPH, Belle and CLEO, compared with the measurement in this paper. The 
statistical and systematic uncertainties are shown separately. 
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Figure 1: The distribution of the mass difference, AM, showing the right-charge 
combinations (points) and wrong-charge combinations (dashed histogram). The 
shaded area shows the signal region, 0.1435 < AM < 0.1475 GeV. The solid line is 
a fit to a Gaussian function plus A{AM — niT^)^ to describe the background, where 
A and B are constants. 
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Figure 2: Distributions of number of D* mesons versus (a) z, (b) p^r^^, (c) 
p^* , (d) rj^* , (e) E^r^^ and (f) rf^^ for data (points) and MC simulations. The 
data are compared with Pythia using the Bowler (solid line) and Peterson, with 
e = 0.06, (dotted line) fragmentation functions and with Herwig (dashed line). 
The component of beauty production as predicted by Pythia (shaded histogram) is 
also shown. 
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Figure 3: Normalised cross section, l/a{da/dz) , for the data (points) compared 
with (a) the symmetric Lund fragmentation modified for heavy quarks (see Eq. 4) 
with vq = 1 (solid line), tq = 0.5 (dashed line) and the original symmetric Lund 
scheme, tq = 0, (dotted line) as implemented in Pythia. The data are also 
compared with (b ) the Peterson fragmentation function with values of the parameter 
e = 0.1 (dashed line), e = 0.01 (dotted line) and the fitted value e = 0.062 (solid 
line) as implemented in Pythia. 
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Figure 4: Normalised cross section, l/a{da/dz) , for the data (points) compared 
with the predictions of FMNRxC^^^ . (a) the Peterson fragmentation function 
in the calculation is shown with e = 0.2 (dotted line), e = 0.035 (dashed line) 
and the fitted value e = O.OTOIqoqq (stat.(Bsyst.) (solid line), (b) the Kartvelishvili 
fragmentation function in the calculation is shown with a = 1.2 (dashed line), 
a = 4.0 (dotted line) and the fitted value a = (stat.®syst.) (solid line). 

The fitted FMNRxC^^^ predictions are shown with the experimental uncertainties 
of the fit (shaded band). 
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Figure 5: D* fragmentation function for the ZEUS data (solid points) compared 
to measurements of the Belle (open circles), CLEO (open triangles) and ALEPH 
(open squares) collaborations in e~^e~ collisions. For shape comparison, the data 
sets were normalised to l/(bin width) for z > 0.3. For the ALEPH data, the 
fragmentation function is measured versus the ratio of the energy of the D* meson 
and the beam energy, whereas for the Belle and CLEO data, the fragmentation 
function is measured versus the ratio of the momentum of the D* meson and the 
maximum attainable momentum at the relevant beam energy. 
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